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This paper presents a developed ecological function for absorption refrigerators with four-temperature- 
level. Moreover, aforementioned absorption refrigerator is optimized by implementing ecological 
function. With the aim of the first and second laws of thermodynamics, an equivalent system is initially 
determined. To reach the addressed goal of this research, three objective functions that the coefficient of 
performance (COP), the ecological function ( E ) and thermoeconomic criterion (F) have been involved in 
optimization process simultaneously. Three objective functions are maximized at the same time. Devel¬ 
oped multi objective evolutionary approaches (MOEAs) on the basis of NSGA-II method is implemented 
throughout this work. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Nowadays, costs of energy consumption and electricity com¬ 
pared to those of fossil fuels have caused engineers in most coun¬ 
tries to consider using absorption refrigerators instead of 
compression refrigerators. However, fuel costs need to be con¬ 
trolled in some way and consequently; economical studies are car¬ 
ried out constantly on optimization of absorption refrigerators and 
reducing their economical costs. 

A new criterion has been used in the recent years to evaluate 
the degree of optimality of thermodynamic cycles. The thermoeco¬ 
nomic criterion economically investigates the phenomena. Capital 
costs and energy costs are considered in these investigations using 
thermodynamic relationships and working conditions are designed 
in a way that the operation is economically optimum. Chen and 
Schouten evaluated the optimal value of coefficient of performance 
in irreversible absorption refrigerator systems [1]. Next, Chen in 
1999 evaluated the optimal value of coefficient of performance 
for the irreversible 4-heat exchange-surface absorption refrigera¬ 
tor in the maximum specific cooling load [2 . This was followed 
by other researches with the aim of optimization of absorption 
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refrigerator processes and other processes similar to those, includ¬ 
ing thermoeconomic optimization of the reversible sterling heat 
pump cycle by Tyagi et al. [3]. The largest portion of capital cost 
of absorption refrigerator pertains to the heat exchangers used in 
the generator, condenser, evaporator, and the absorber. 

From economical point of view, incorporation of absorption 
refrigerators with the combined power and heat (CHP) systems 
has high profit. Since the 1970s finite-time thermodynamics have 
been used for thermal systems analysis. Angulo-Brown was sug¬ 
gested an ecological criterion [4] for optimizing refrigerator sys¬ 
tem. The aforementioned suggested criterion is illustrated as 
E = W - T l g for the finite-time Carnot machines. In which, T L 
stands for the temperature of cold reservoir, and a represents the 
production of entropy production throughout the system. Yan ver¬ 
ified which the ecological index as E = W - T 0 & is more practical 
[5]. Where, T 0 is the ambient temperature. There are a great num¬ 
ber of attentions in the previous works for the optimization of irre¬ 
versible absorption refrigerators. Bhardway and colleagues 
optimized absorption refrigerators utilizing finite-time thermody¬ 
namics [6,7 . Moreover, Huang and colleagues studied the COP, 
surface area, and ecological principles for the irreversible absorp¬ 
tion with four-temperature-levels [8]. Chen presented models for 
the optimization of absorption refrigerators [9,10]. Sun and col¬ 
leagues investigated the optimization of the relationship between 
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Fig. 1 . Generalized four-temperature-absorption refrigerator system [3]. 


entropy production and heating load of the internal heating heat 
exchanger [11]. Furthermore, Chen and colleagues presented an 
approach for the analysis of absorption refrigerator cycles with 
four-temperature-levels [12 . Chen and Yan optimized the ecolog¬ 
ical criterion of an irreversible absorption heat exchanger [13]. 

Two well-known researchers [14,15] optimized the absorption 
cycles for the ecological criterion for the irreversible Carnot cycles 
with three-temperature-levels. Fathi et al. presented an approach 
for solar and irreversible absorption cooling [16]. Furthermore, 
excluding for the ones specified now, there are numerous optimiza¬ 
tion of performance investigates for absorption refrigerators 
[17-20]. Moreover, latest optimization approaches named 
‘advanced exergy based analysis approaches’ [21,22] or ‘exergy- 
entropy relationship’ [23] of the system can be studied for the 
further researches. Hellmann [24] analyzed the effect of a supple¬ 
mentary restriction, or of an additional operational variable 
required. Flajtaleb and Feidt, and Vasilescu et al., have also evolved 
sensitivity analysis to a 4 heat reservoirs refrigerator [25,26]; these 
two communications were completed [26] by optimization of the 
COP for a four finite reservoirs reverse cycle irreversible machine 
(heat pump or refrigerator), and the analogous parametric analysis 
was accomplished. 

Multi-objective optimization is useful way for solving engineer¬ 
ing problems at various fields for research [27-29 . 

Solving the optimization of multi-objective problem requires 
synchronized satisfaction of different objectives. Thus, evolution¬ 
ary algorithms (EA) were established and evolved through the 
20th century in an effort to solve multi-objective issues at random 
[30]. An suitable way to a multi-objective quandary is to inquiry an 
collection of routes, each of them convinces the objectives at a sat¬ 
isfactory degree away being overcome by another route [31]. 
Multi-objective optimization quandaries generally stand for a fea¬ 
sibly countless group of routes which named frontier of Pareto, 


whose assessed vectors show the foremost feasible exchanges 
throughout the objective function region. Recently, multi-objective 
optimization of diverse thermodynamic and energy systems were 
employed in different field of science [32-45 . 

In our study, inorder to maximize the coefficient of performance 
(COP) and the ecological function ( E ) and thermoeconomic crite¬ 
rion (F), multi-objective optimization algorithms has been 
employed. Error analysis was performed to determine accuracy 
of ultimate results of various decision making methods. 

2. Thermodynamic analysis of system 

Thermodynamic analyses have been performed for the equipol¬ 
lent system, and earlier equations presented by Refs. [46-53] have 
been implemented throughout the determinations. The ambient 
temperature was supposed as 290 K. The scheme of the general 
absorption refrigerator which under studied here is depicted in 
Fig. 1. 

2 A. Thermodynamic analysis of the generator-absorber assembly 
(heat engine part) 

Thanks to the first law of thermodynamics, 

Qg-(Ia-W = 0 ( 1 ) 

In which, Qa (heat output to the absorber heat sink) and Q c 
(heat input from the generator heat source) can be formulated as 
following as: 

0.G = ocA g (Tc — T g ) (2) 

Q_ A = PAA(T a -T A ) (3) 

Thanks to the second law of thermodynamics [13-15,46-53], 
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In which, I H represents the parameter of internal irreversiblity 
for generator-absorber assembly. I H = 1 should be consider for 
the reversible cycles, and I H > 1 for the irreversible cycles. Heat 
transfer surface area ratio (a H ), total heat transfer surface area 
(. A t ,H ) and Working fluid temperature ratio (x H ) for the absorber- 
generator assembly can be defined as following [53]: 


x H 




a H 


Ai 


A ( 


( 6 ) 


Atm — Ac + A& 


(7) 


T g , T a , and Q.g can be formulated by employing Eqs. (2)-(7) as 
following as [53]: 


T — 

1 g — 


T a = 


Qc = 


I h ocT g x h + cihPTa 

x h {cihP + Ih&) 

(8) 

IhOcTcXh + OhPTa 
( a H f + I H oc) 

(9) 

%At,h (a H P(T G x H — Ta)\ 

1 + % V x H (a H p + I H oc) J 

(10) 


2.2. Thermodynamic analysis of the evaporator-condenser assembly 
(refrigerator part) 


Based on the concept of the first law of thermodynamics, 

Qg-Qe- W = 0 (11) 

Here, Q.e (heat input to the evaporator) and Qc (heat rejection from 
the condenser) are 

Qc = JuA c (T c - T c ) (12) 

Qe = 9A e (T e - T e ) (13) 

According to the second law of thermodynamics [13,15,46-53], 


Qc 

T c 

Qc 

T c 




In which, I R represents the parameter corresponds to internal 
irreversiblity for evaporator-condenser assembly. I R = 1 should 
be consider for the reversible cycles, and I R > 1 for the irreversible 
cycles. Total heat transfer surface area (A TtR ), and x ratio of heat 
transfer surface areas (a R ), and total heat transfer surface area 
( A T}R ) for the evaporator-condenser assembly can be defined as fol¬ 
lowing formulas [53]: 


T 

1 e 

X R =y 

1 c 

(15) 

Ac 

aR ~T r 

(16) 


At,r — Ac + Ae 

T c , T e , Q.c and Qe can be defined using Eqs. (12)-(17): 
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j _ pT G x R a R + 6I r T r 
c ~ x R (a R [i + I R 0) 

j _ fiT G x R a R + 9I r Te 
e ~ (a R ii + I R 0) 

■ _ 6Aj R n R a R p(TE — TcX r ) \ 

1 + V x R (a R p + l R 0) J 


(18) 

(19) 

( 20 ) 


• _ 6At R f a R p(T E — TcX r ) \ 

1 + a R \ ( CL R p + I R 0) J 


( 21 ) 


The evolved ecological function which includes the magnitudes 
of maximum energy destruction of the system, maximum power 
output for the heat engine, and maximum cooling load in the 
absorption cooling system is defined as following as [46-53]: 


r n /"\ o nr I 0-A Qg Qc Qf 

T Qe Qg ~ Qa ~ T 0 ( -^ t- ~t~ 

IA * G 1 C i E 


( 22 ) 


The coefficient of performance of the absorption system can be 
formulated as following as [13,15,49-53,55]: 


COP = 


Qe 

Qg 


(23) 


And the specific cooling load is obtained as follows: 

Qr 


A 


(24) 


A — Atm + (24a) 

Eqs. (23) and (24) are the general relationships for absorption 
refrigerators with four reversible heat sources. These relationships 
can be used to obtain thermo-economical optimum performance 
for this type of absorption refrigerators. According to the studies 
carried out by Kodal et al. as given in [54 , the thermoeconomic cri¬ 
terion for absorption refrigerators with four reversible heat sources 
is defined as the total price of unit cooling load which includes 
both capital and energy costs. Therefore, the function which is to 
be optimized is as follows: 




where C„ C e are capital and energy costs in the unit of time, respec¬ 
tively. The capital cost of absorption refrigerators is assumed to be 
proportional to its total heat transfer area: 

Ci = k\A (26) 


where /<t is the capital recovery factor, the capital cost for the unit 
heat transfer area. Energy consumption cost is directly proportional 
to the rate of heat input: 



= foO.G 



where k 2 is the unit cost of energy. By substituting Eqs. (26) and 
(27) in Eq. (25), we obtain: 


/ 


Qe 

/<iA + /<2Qg 


1 

hi _i_ Ti _ 
r ^ COP 



By defining k = ^ as the thermoeconomic parameter having the 
dimension when the capital cost increases and the cost of 
energy consumption reduces, the thermoeconomic parameter k 
goes up. 

The ration of the optimal value of the thermoeconomic criterion 
of absorption refrigerators with four reversible heat sources is 
obtained as follows: 
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Fig. 2. Pareto frontier (Pareto optimal solutions) using NSGA-II. 
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Fig. 3. The distribution of a H for the optimal points on Pareto front. 
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Fig. 4. The distribution of a R for the optimal points on Pareto front. 

(29) 


3. Optimization process 

3.1. Evolutionary algorithm: genetic algorithm 

Genetic algorithms use a repetitive, random search method to 
find the best answer and copy it in an easy way principle of biolog¬ 


ical evolution [56 . First parameter in an evolutionary algorithm 
would be the population of individuals. Individual expresses the 
values of the decision variables and is a potential answer to the 
optimization problem [56 . More about genetic algorithm and its 
process could be finding in references: [56,57]. 


3.2. Objective functions, decision variables and constraints 

The coefficient of performance (COP) and the ecological function 
(E) and thermoeconomic criterion (F) three objective functions for 
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Population 


Fig. 5. The distribution of x H for the optimal points on Pareto front. 



Fig. 6. The distribution of x R for the optimal points on Pareto front. 


Table 1 

Decision making of multi-objective optimal solutions. 


Decision making method 

Decision variables 



Objective functions 



dR 


Xr 

COP 

E (kW) 

F 

TOPSIS 

2 

0.3 

0.85 

0.590 

4.33 

908.6 

2.75 

LINMAP 

2 

0.3 

0.85 

0.587 

4.54 

829.4 

2.89 

Bellman Zadeh 

2 

0.3 

0.85 

0.565 

4.65 

784.2 

2.96 


this work, which are identified by Eqs. 22, 23 and 29. Four decision 
variables have been choosed in our study as follows: 

• a H : ratio of heat transfer surface area. 

• a R : ratio of heat transfer surface area. 

• x H : hot working fluid temperature ratio. 

• x R : refrigerant working fluid temperature ratio. 

Albeit the decision parameters would be different in the optimi- 


zation process, each should be in an appropriate interval. By con¬ 
sidering following limitations the objective functions was 
unraveled: 

2 ^ Qh ^ 2.5 

(30) 

0.2 ^ a R ^ 0.3 

(31) 

0.85 ^ Xj-i ^ 0.95 

(32) 

0.55 ^ x R ^ 0.6 

(33) 


4. Results and discussions 

Utilizing multi-objective optimization on the basis of the 
NSGA-II method, both the coefficient of performance (COP) and 
the ecological function (E) and thermoeconomic criterion (F) are 
maximized at the same time. The objective functions in the per¬ 
formed optimization, and the limitations which have been used, 
are identified by Eqs. 22, 23 and 29 and Eqs. (30)-(33) respectively. 
The Working fluid temperature ratio (x H , x R ) and heat transfer sur¬ 
face area ratio (a H , a R ) are assumed as design parameters through¬ 
out the optimization process. Due to agreement with appreciated 
literatures, following characteristics of the refrigerator system are 
presumed as below [46-53]: 

A e = 1 50, A c = 45, I H = 1.05, I R = 1.05, T c = 313 (K), T E = 288 (I<), 

T g = 393 (K), T 0 = 290 (K), a = 2.3, p = 0.4, p = 3.265, 9 = 0.195, T A = 
315 (K). 

Fig. 2 demonstrates the frontier of optimum Pareto for objective 
functions E and COP and F, also optimal outcomes of TOPSIS, 
LINMAP and FUZZY decision making techniques. The span of 
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Fig. 7. Effects of the ratio of heat transfer surface area ( a R ) and the ratio of heat 
transfer surface area (a H ) on the (a) coefficient of performance, (b) ecological 
function, and (c) thermoeconomic criterion. 


Fig. 8. Effects of the ratio of heat transfer surface area ( a R ) and the hot working fluid 
temperature ratio (x H ) on the (a) coefficient of performance, (b) ecological function, 
and (c) thermoeconomic criterion. 
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Fig. 9. Effects of the refrigerant working fluid temperature ratio (. x R ) and the hot 
working fluid temperature ratio (x H ) on the (a) coefficient of performance, (b) 
ecological function, and (c) thermoeconomic criterion. 


Fig. 10. Effects of the refrigerant working fluid temperature ratio ( x R ) and the ratio 
of heat transfer surface area ( a H ) on the (a) coefficient of performance, (b) ecological 
function, and (c) thermoeconomic criterion. 
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Table 2 

Error analysis based on the mean absolute percent error (MAPE) method. 


Error (%) 

TOPSIS 



LINMAP 



Fuzzy 




COP 

E 

F 

COP 

E 

F 

COP 

E 

F 

Max error 

0.4 

1.1 

0.4 

0.9 

2.0 

0.9 

1.7 

2.6 

1.7 

Average error 

0.2 

0.6 

0.2 

0.2 

0.5 

0.2 

0.7 

1.1 

0.7 


variations for optimum outcomes of F is from 682.6 to 1057 (kW), 
and from 3.543 to 4.848 for COP and F is from 2.254 to 3.085. 

Figs. 3-6 depict the variation of different magnitudes of deci¬ 
sion parameters throughout their accepted span for the optimum 
design points on the front of Pareto. The variation of the a H variable 
is illustrated in Fig. 3. 

In Fig. 4 the blue line, shows the greatest magnitude of a R which 
is equal to 0.3. Fig. 5 shows that, the variation of different magni¬ 
tudes of x H for the optimum design points on the front of Pareto is 
equal to 0.85. The variation of x R for the optimum points on the 
front Pareto is depicted in Fig. 6. 

Table 1 shows the optimum outcomes gained for objective 
functions and decision variables employing Fuzzy Bellman-Zadeh, 
TOPSIS and LINMAP techniques. 

4.1. Sensitivity analysis 

From Fig. 7 it can be seen that, the coefficient of performance 
and ecological function and thermoeconomic criterion increase 
considerably with increasing of the ratio of heat transfer surface 
area (a R ) at various values of the ratio of heat transfer surface area 

In Fig. 8 different system parameters against the ratio of heat 
transfer surface area ( a R ) in the different hot working fluid temper¬ 
ature ratio ( x H ) are plotted. And in these figures by increasing the 
ratio of heat transfer surface area ( a R ), the coefficient of perfor¬ 
mance and ecological function and thermoeconomic criterion will 
be raised. 

In Fig. 9 different system parameters against the refrigerant 
working fluid temperature ratio ( x R ) in the different hot working 
fluid temperature ratio ( x H ) are shown. As can be seen, with 
increasing the refrigerant working fluid temperature ratio, the 
coefficient of performance and thermoeconomic criterion decrease 
and ecological function will be raised. 

In Fig. 10 different system parameters against the refrigerant 
working fluid temperature ratio (x R ) in the ratio of heat transfer 
surface area ( a H ) are shown. As can be seen, with increasing the 
refrigerant working fluid temperature ratio ( x R ), the coefficient of 
performance and thermoeconomic criterion decrease and ecologi¬ 
cal function will be raised. 

Finally, to determine deviation of the employed methods, error 
analysis is used. Error analysis is performed according to average 
absolute relative deviation (AARD) and maximum absolute relative 
deviation (MARD) of the outcomes obtained from each decision 
making methods. It should be noted that, error analysis is applied 
on the final solutions which gained after 30 runs of aforemen¬ 
tioned decision making methods. Table 2 summarized outcomes 
of error analysis for each decision making approach. 


5. Conclusions 

With the aim of evolved ecological function, the optimization of 
irreversible absorption refrigerator with four-temperature-levels 
has been carried out. A balance between entropy production and 
the load supplied to the system and cooling load has been estab¬ 


lished. The effects of x H , a R , x R , and a H factors and temperature of 
operation fluids on the aforementioned system are analyzed. 

The ecological function (F), the coefficient of performance (COP) 
and thermoeconomic criterion (F) are simultaneously assumed for 
multi-objective optimization, in order to determination by ther¬ 
modynamic analysis. The ratio of working fluid temperatures 
( x H , x R ) and ratio of heat transfer surface areas ( a H , a R ) have been 
presumed as design variables. The multi objective evolutionary 
algorithm on the basis of the NSGA-II method has been utilized 
and the frontier of optimum Pareto throughout the objectives 
space, was specified. To specify an ultimate solution from the 
results gained by multi-objective optimization process, three 
well-known decision making techniques have been implemented 
and final outcomes are compared on the basis of error analysis. 
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